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Abstract Using a suite of simulations Q which sucessfully produce bulgeless 
(dwarf) disk galaxies, we provide an analysis of their associated cold interstellar 
^^ i media (ISM) and stellar chemical abundance patterns. A preliminary comparison 

rj I with observations is undertaken, in order to assess whether the properties of the 

cold gas and chemistry of the stellar components are recovered successfully. To this 
end, we have extracted the radial and vertical gas density profiles, neutral hydrogen 
velocity dispersion, and the power spectrum of structure within the ISM. We com- 
2 I plement this analysis of the cold gas with a brief examination of the simulations' 

metallicity distribution functions and the distribution of a -elements-to-iron. 
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rr' ' 1 Introduction 

QQ ■ Historically, simulating the formation and evolution of a bulgeless, rotating, stel- 

^f) ', lar disk, within the context of the classical picture of hierarchical assembly, has 

VO ' proven problematic [1]. Overcoming this has only been achieved recently [2J, with 

\l^ I Governato et al. presenting simulated dwarfs with bulge-to-disk ratios <0.04, with 

f"~-« ■ very shallow associated central dark matter profiles. This was achieved by resolv- 

^D I ing the inhomogeneous ISM, thereby allowing for the imposition of a realistic star 

formation density threshold which resulted in strong supernovoae-driven outflows 
preferentially removing low angular momentum gas. 

In what follows, we concentrate our analysis on three simulated dwarf galaxies, 
the basic properties for two of which (DGl, DGILT) have been introduced pre- 
viously |I2- In fj2l we quantify the spatio-kinetic characteristics of these bulgeless 
5^ \ dwarfs, and contrast them with recent observational work |l4l|5]|6]. A preliminary 

analysis of the chemical properties of the associated stellar components is then pre- 
sented in fj3]and compared with relevant empirical work in the field QUI]. As the 
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simulations are not modelled on any one particular galaxy, the analysis should be 
seen as comparative only in a qualitative sense. 



1.1 Simulations 

The three simulations analysed were run using the N-body+SPH code GASOLINE 
ifTOll . employing the same initial conditions, and differing only in the adopted res- 
olution, star formation density threshold, feedback efficiency, and star formation 
efficiency. The base simulation (DGl) III uses a star formation density threshold a 
thousand times higher (100 cm^^) than that adopted in previous cosmological simu- 
lations (~0.1 cm^^); a lower density threshold version (DGILT), despite its posses- 
sion of classical problematic traits (i.e., centrally-concentrated light (B/D«0.3) and 
dark matter density distributions), is included here for comparison. This implemen- 
tation of a more realistic density threshold for star formation was made possible by 
the associated high force resolution (86pc) which allowed us to resolve individual 
star forming regions of mass ~'10^ M^,. We include in our analysis a new simula- 
tion (nDGl), similar to DGl, but with two important modifications - the inclusion of 
high-temperature (T>10^K) metal-line cooling [9], and enhanced supernova feed- 
back (100% of the energy being coupled to the surrounding ISM in the form of 
thermal energy, rather than the 40% used in the base simulation). The feedback 
method is described fully in [8J . 



2 HI Analysis of the Simulated Dwarfs 

Throughout our analysis, we cross-check results which are based on cold gas parti- 
cles (T<15000K) i8| with those based on the inferred neutral hydrogen, to ensure 
our comparisons with empirical data have not been biased in any obvious manner. 
As the current implementation of ISM physics within GASOLINE does not allow 
gas to cool to densities more appropriate to molecular hydrogen, what we label as 
'HI' can possess column densities '^S-lOx larger than encountered in nature. For 
that reason, one must be careful to not overinterpret the quoted surface densities. 



2.1 Radial Density Profiles 

We first examine the radial density profiles of the cold gas in the three simulations 
(Fig [Hi- Observations of HI disks 161 14| show that, in nature, they extend from ^2 
to ^6 radial disk scalelengths and truncate near ^1.5r^. DGILT shows a flat profile 
with a formal scalelength of ~18 kpc, which truncates at ^O.Sr^/. The base simula- 
tion (DGl) has a radial scalelength of ^6 kpc (which is also approximately where 
it truncates); embedded within this extended disk is a high density cold gas core. 
Our new simulation (nDGl) can really only be fit using two exponentials, spanning 
the inner and outer disks (such double exponentials are fairly common in nature), 
and also possesses a compact cold core. These cores are somewhat transient, being 
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disrupted during periods of more significant star formation; e.g., for DGl, at z^Q-7 
(the last epoch of substantial star formation - O.lMQ/yr compared to 0.005Mp/yr at 
the present-day), the cold core is absent. These conclusions are robust to the choice 
of 'HI', as opposed to 'cold gas' . 



O DGl Cold Gas (r^ = 5.7kpc) 
XnDGI Cold Cos (r^i = 2.3kpo) 
X nDGI Cold Gas (rj'„ = Q.9kpc) 
A DGl LT Gold Gas (r^' = 1 8.3kpc) 




Radius [kpc] 



Fig. 1 Radial cold gas density profiles for the simulated dwaifs DGl (diamonds), DGILT (trian- 
gles), and nDGI (crosses). The thick overplotted lines show the exponential fits to the distributions, 
from which the noted scalelengths were derived. nDGI is best represented by a double exponential, 
with a break between the two near ~3 kpc. The cold gas of DGl is distributed in a more extended 
exponential disk component of scalelength ~6 kpc, while that of DGILT is ~18 kpc. 



2.2 Flaring 

Both DGl and nDGI, when viewed edge-on, show significant flaring of their gas 
disks. In nature, comparable dwarfs |4| typically show an increase in the FWHM 
of their vertical density distribution of ^50% in going radially from QAid to 1 .Or^. 
Fig ID shows the vertical density profiles for the simulations at three annuli (O.h^, 
QAiii, and l.Or,/); we can see that DGl flares by a factor of ^4 from 0.4rf/ to l.Or^, 
while nDGI flares by a factor of ~1.5. In that 'fractional' sense, DGl seems some- 
what extreme, but quantifying the flaring in terms of physical units (kpc), the degree 
of flaring is not dissimilar to that observed 1|4J. 



2.3 Velocity Dispersion 

Observations of HI velocity dispersion profiles |6|, for both giants and dwarfs, 
show that galaxies possess characteristic dispersions of ^8 — lOkm/s (rising to 
'^12—15 km/s in the inner star forming part of the disk). In Fig [3] we show the 
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Fig. 2 Vertical density pro- 
files of tlie cold gas disks 
for nDGl (upper panel) and 
DGl (lower panel), where 
'cold' here means gas parti- 
cles with temperatures less 
than 15000 K. The profiles 
are measured in three seperate 
annuli corresponding to 0.1, 
0.4, and 1.0 HI disk scale 
lengths. Overlaid on the plots 
are the exponential fits from 
which we derived the asso- 
ciated scaleheights (inset to 
each panel). 
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HI line-of-sight velocity dispersion profiles for nDGl, DGILT, and Holmberg II 
161. Critically, we show nDGl with (squares) and without (crosses) the inclusion of 
thermal broadening, following the methodology of I.12J . Because SPH tracks only 
the streaming motions of the gas particles, one must incorporate the random veloc- 
ity component to each particle, using the internal energy of each particle as input. In 
practice, we draw random velocities for each Cartesian coordinate from a Gaussian 
with CT = \JkT j \x and add those to the components of the streaming motion. With- 
out the inclusion of thermal broadening, nDGl (and DGl) shows an extremely (and 
unphysically) kinematically cold ISM compared to that of DGILT (which despite its 
aforementioned problems, resembles empirical data such as HoII very well). Two 
unresolved issues which are apparent from Fig |3] are that (i) while the dispersion 
based upon the streaming velocities (crosses) shows the same increase in the star 
forming portion of the disk that is observed in nature ||6l, the magnitude of the ther- 
mal broadening essentially "wipes out" this signal, and (ii) the velocity ellipsoid 
of the cold gas becomes isotropic, disguising any anisotropics that might have been 
present in the streaming motions (i.e., young stars, and the cold gas from which they 
formed, will necessarily have different velocity ellipsoids). 

When DGl is viewed at redshift z~0.7, the more turbulent ISM being driven by 
the enhanced and sustained star formation (0.1M©/yr compared to O.OOSo/yr), was 
a close match to the velocity dispersion profiles observed in nature - put another 



Chemodynamical Simulations of Dwaif Galaxies 5 

way, the dispersion of the streaming motions during periods of significant star for- 
mation are already on the order of the dispersions observed in nature, without the 
application of thermal broadening. As the gas clouds out of which stars are form- 
ing within these simulations have temperatures on the order of '--^SOOOK, there is 
no way to avoid the imposition of a significant ^8 km/s thermal broadening term; 
until we can resolve densities corresponding to cores of molecular clouds, this SPH 
limitation must be observed. 



2.4 Power Spectrum 

We next derived the spatial power spectra of the HI gas associated with our simula- 
tions, as well as that for the Small Magellanic Cloud |5|, as shown in Fig 21 using 
Fourier Transforms of the moment zero density maps, after |5|. In the simulated 
galaxies, we smoothed with a lOOpc Gaussian to be consistent with the beam smear- 
ing present in the observational data. All of the galaxies were then fit with a power 
law PocF, where 7=-3.5 for DGl, 7=-3.4 for DGILT, 7=-4.2 for nDGl, and 
7=— 3.2 for the SMC. The steeper slope possessed by nDGl indicates an excess of 
power distributed on large scales within its simulated ISM, relative to that observed 
in the other simulations or the SMC (i.e., the enhanced feedback has shifted power 
from smaller scales to larger ones). The observational data (for the SMC and other 
dwarfs) is well fit by a pure power law, indicative of a lack of preference in HI cloud 
size, in nature), while the simulations show departures from this scenario. These 
departures can be traced to specific structures in the simulations - e.g., in nDGl, the 
enhanced power seen on scales of ^^400— 500pc corresponds to the spacing of the 
tightly-wound spiral structures seen in its inner few kpcs. 



3 Chemistry of the Simulated Dwarfs 

The next direction for our work will be based upon a comparison of the stellar chem- 
istry of the simulations with existing observational data. To date, we have exam- 
ined the metallicity distribution functions (MDFs) and oxygen-to-iron distributions. 
Here, we simply highlight preliminary results pertaining to nDGl. 

Fig|5]shows the present-day cumulative [Fe/H] MDF for all stars within 4 kpc of 
the three simulations, normalised to unity (arbitrarily) at [Fe/H]=— 2.3. Focusing on 
the metal-poor tail, it is apparent that our simulated dwarfs suffer from a "G-dwarf 
problem" (relative overproduction of metal-poor stars), an issue to which we will 
return in an upcoming study. 

The global MDFs of the simulations are similar to those seen in Local Group 
dwarfs, with an ^0.16 dex dispersion 'peak' at [Fe/H] s»— 1.1 superimposed on a 
broader metal -poor component (the formal dispersion of entire nDGl MDF be- 
ing ~0.32 dex). This is not dissimilar to the ^^0.25 dex (intrinsic) dispersion seen 
(for example) in the Local Group dwarf. Carina fS\. While nDGl is not a dwarf 
spheroidal, it does have an episodic star formation history akin to that of Carina. 
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Fig. 3 HI line-of-sight (40° viewing angle) velocity dispersion profiles of nDGI (crosses: particle 
streaming motions only), nDGI with the inclusion of thermal broadening (squares), DGILT (tri- 
angles), and observational data for Holmberg II (plus signs), scaled to ^25 in the B-band (isophotal 
radius corresponding to 25 mag/arcsec", roughly the extent of the star forming disk). 
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Fig. 4 Power spectra of the HI distributions of the three simulations - nDGI (crosses), DGl (di- 
amonds), and DGILT (triangles) - accompanied by that for the SMC. We truncated the power 
spectra for the simulations at ~200pc (twice the beam size used to simulate the effects of beam 
smearing). The break in the SMC spectrum is due to a missing base line fSJ. Overplotted on each 
dataset is a representative power law fit to the data, with the slopes noted in the inset to the panel. 
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Fig. 5 Cumulative MDFs for nDG 1 (dotted line), DO 1 (dashed line) and DO ILT (dot dashed line). 
The compilation of observational data shown 1 7 1 includes that of the Galactic halo (asterisks) and 
the combined MDF constmcted by merging those of the Carina, Sextans, and Sculptor dwarfs from 
DART data using two different metallicity calibrations (triangles and diamonds). 



Fig|6]shows the relationship between [O/Fe] and [Fe/H] for the same stars (black 
dots) within 4 kpc of nDGl. Overplotted, again not because it is supposed to be a 
true analog of nDGl but because it provides a useful benchmark, are the [a/Fe]- 
[Fe/H] data (asterisks) for the Local Group dwarf. Sculptor ifTTl . The distribution 
of the residuals in [O/Fe] about the best-fit lines through the Sculptor and nDGl 
datasets are both consistent with intrinsic scatters of ^0.13 dex. For nDGl, this 
scatter varies somewhat with metallicity, with the scatter in the [O/Fe] residuals 
for stars near [Fe/H]sa— 1 being ^0.1 dex, while those with [Fe/H] < — 1.5 show a 
scatter closer to ^0.2 dex. To first order though, it would suggest that the adopted 
magnitude of metal diffusion employed was reasonable 
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Fig. 6 [0/Fe]-[Fe/H] distribution for all stars within 4kpc of nDGl. For comparison, data from 
the Local Group dwarf, Sculptor ! ifl , are overplotted (asterisks). 



